Noise-induced effects in high-speed reversal of magnetic dipole 
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The effect of noise on the reversal of a magnetic dipole is investigated on the basis of computer 
simulation of the Landau-Lifshits equation. It is demonstrated that at the reversal by the pulse with 
sinusoidal shape, there exists the optimal duration, which minimizes the mean reversal time (MRT) 
and the standard deviation (jitter). Both the MRT and the jitter significantly depend on the angle 
between the reversal magnetic field and the anisotropy axis. At the optimal angle the MRT can be 
decreased by 7 times for damping a=l and up to 2 orders of magnitude for a=0.01, and the jitter 
can be decreased from 1 to 3 orders of magnitude in comparison with the uniaxial symmetry case. 



The decrease of sizes of magnetic nanoparticles used in 
storage media leads to the increase of fluctuations and, 
therefore, to increase of storage and switching errors (jit- 
ter). Thus, theoretical investigation of noise-assisted 
high-speed switching of magnetic dipoles is of crucial 
importance. The most studies are based on computer 
simulation of the Landau-Lifshits equation^ describing 
the dynamics of the magnetic dipole in a magnetic field. 
It has been found that there is an optimal angle between 
the applied magnetic field and the anisotropy axis, which 
is typically around 45 degrees^^; with the decrease of 
the magnetic field rise time, the coercivity of magnetic 
particle (dynamic coercivity) increases^,—. The influence 
of the damping and the external magnetic field ampli- 
tude on the reversal time has been investigated \v^^^; 
the dependence on the pulse shape of the field has been 
studied in^. All the above mentioned results have been 
obtained at zero temperature, i.e. without account of 
thermal fluctuations. However, at finite temperature the 
most important problem is the stable magnetic reversal: 
the remagnetization process must occur with minimal 
switching time and the standard deviation. Unfortu- 
nately, there were very little number of works, devoted 
to investigation of statistical characteristics of magnetic 
reversal processes. The mean reversal time (MRT) has 
been studied in Ref.sS.,-", and in Ref.^° it has been shown 
that during noise-assisted reversal, the noise leads to the 
decrease of MRT. 

In the present paper the investigation of the rever- 
sal process of a magnetic dipole has been performed on 
the basis of computer simulation of the Landau-Lifshits 
equation with thermal fluctuations taken into account. 
It is focused on the investigation of statistical charac- 
teristics of the reversal process with the aim to find an 
optimal regime of reversal with the smallest mean rever- 
sal time and the standard deviation. 

The dynamics of magnetic dipole is described by the 
Landau-Lifshits equation: 
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where a} is the magnetization of a particle, 7} is the 
effective magnetic field, 7 is the gyromagnetic constant. 



is the saturation 



f3 = a is the damping, Mg 

magnetization. The effective magnetic field contains the 
following components: ^ — + hI + Hj-, where 11^ 
is the anisotropy field. He is the external field, and Ht 
- fiuctuational field. The fluctuational field is assumed 



to be white Gaussian noise with zero mean and the cor- 
relation function: {H{t)nH{t')Tj} = ^w^S{t-t')6ij, 
where k - Boltzmann constant, T is the temperature, 
and V is the volume of the magnetic particle. 

Let us consider the reversal of magnetic dipole, ini- 
tially magnetized along anisotropy axis and along x-axis 
from the state a}[+Ms,0,0] to the state A?[-Ms,0, 0]. 
To find the area of parameters where the fastest and the 
most reliable reversal occurs, as the characteristic to be 
studied let us choose the first passage time of a certain 
boundary. The mean first passage time (the mean rever- 
sal time, MRT) t, and the standard deviation of the first 



passage time cr (SD, jitter) are 



11. 



where t 



is the first 



passage time of an absorbing boundary and N > 10000 
is the number of realizations. As in Ref.^ let us choose 
the boundary as the passage of the point A?[0, My, M^]. 

In the calculations it is convenient to use the param- 
eters, related to the magnetic recording media^^: a=0.1, 
7=1.76 • lO^Hz/Oe, M^=360 Oe, V=2 ■ lO^nm^, the 
anisotropy constant K—7.2 ■ lO^erg/cm'^. The static co- 
ercivity is Hc—2K/Ms='lQ00 Oe. For modeling we take 
the amplitude of the magnetic field to be _ffo=6000 Oe. 
It is known that the driving by the signal with sharp 
fronts leads to the minimal MST^^. However, the pulses 
used in real recording media systems have finite rise 
time^. As an example of a driving with smooth fronts 
we consider the sinusoidal pulse He=eHQshnTt/tp with 
the width tp (see the inset of Fig. 1), where e is unitary 
vector of magnetic field direction. If the switching dur- 
ing tp does not happen, the computation is continued 

for He=0 until some long period of time tf, much larger 
than any other relaxation time scale. Our aim is to find 
the parameter range, for which the reversal by a smooth 
pulse takes place quickly and reliably. 

The Landau-Lifshits equation with noise has been 
computed both by the Heun method programmed in 
Fortran and by the specialized package SIMMAG (Sim- 
ulation of MicroMAGnets), developed in the laboratory 
of mathematical modeling of IPM RAS. 

It is known that for zero temperature T — the re- 
versal of the dipole by the longitudinal field, 6 = 0, does 
not occur, since the dipole is in the equilibrium state, 
even if this state is unstable. The presence of thermal 
fluctuations allows to move the dipole away from this 
equilibrium state. In Fig. 1 the plots of the mean rever- 
sal time T and the standard deviation a are presented. 
First of all, it is seen that both r and cr have minima as 
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FIG. 1: The mean reversal time (solid curves) and the stan- 
dard deviation (dashed curves) versus pulse width for the 
zero angle between anisotropy axis and the driving field 9=Q. 
Inset: the driving pulse. 
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FIG. 2: The MRT (solid curves) and SD (dashed curves) 
versus pulse width for different angles between anisotropy 
axis and external field and different temperatures for 6 = 5°. 



functions of the driving pulse width. This indicates that 
both these temporal characteristics can be minimized by 
the optimal choice of pulse duration. Similar effect has 
recently been observed for Josephson junctionsi^. The 
decrease of the MRT at large durations is due to the fact 
that with decrease of the width the potential barrier dis- 
appears faster. With further shortening of the pulse, the 
magnetization does not have enough time for the com- 
plete reversal during tp, so the MRT increases. This, 
actually, means that for rather short pulses the transi- 
tion occurs due to effect of fluctuations (the so-called 
noise-induced switching). The standard deviation with 
decrease of the pulse width behaves similarly to r, but 
the curves for different temperatures cross each other. 
For long pulses smaller temperature leads to smaller a, 
but for short pulses - vice versa. This is also explained 
by the transition from the regime of switching by the ex- 
ternal field to the noise-induced escapes, since it is well- 
known, see, e.g.jii, that in the noise- induced regime the 
SD is approximately equal to the MRT. 

Let us point out few more peculiarities, which are 
clearly visible in Fig. 1: higher temperature leads to 
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FIG. 3: The MRT (solid curves) and SD (dashed curves) ver- 
sus the angle between anisotropy axis and external magnetic 
field for different values of damping a, tp = 8ns and temper- 
ature T=300 K. Inset: trajectories of the magnetization for 
different angles, a=0.1, tp=15.7 ns, T=300 K; 9 ^ 0° - blue 
curve, d = 5° - red curve, d = 45° - green curve. 
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FIG. 4: The MRT (sohd curves) and SD (dashed curves) ver- 
sus the angle between anisotropy axis and external magnetic 
field for different values of magnetic field amplitude, tp=8 ns 
and temperature r=300 K. 



smaller MRT, i.e. noise allows to speed up the reversal, 
which agrees with the predictions of ReffiS, but contra- 
dicts with the results of Ref.-l^. Besides, the dependence 
a ^ VT usual for Josephson junctions at large pulse 
widths, is not also reproduced here, from Fig. 1 one 
can see that (t{T) dependence is slower than \/T- The 
explanation of such an unusual behavior is the specific 
nature of the reversal process at 9 — 0: in difference 
with the Josephson junctions^^, due to location of initial 
condition at the unstable equilibrium point, the reversal 
is impossible at zero temperature and the deterministic 
reversal trajectory does not exist. That is why fluctua- 
tions help to leave the unstable initial state and namely 
this leads to the described above peculiarities. It should 
be noted that the above investigated case 9^0, which 
is mostly studied analyticall y^"^ is degenerate: first, 
technically the angle between the anisotropy axis and 
the external field can be set up with a certain precision; 
second, the MRT is largest in this case, see below. 



In Fig. 2 the MRT and SD are given for three different 
values of angles 0=0°, 5°, 45° for T=300 K and three 
different temperatures for 9=5°. First, let us focus on 
the curves for 9=5°. In spite that at large tp one can see 
little decrease of MRT with increase of the temperature, 
at small tp around minimum the opposite effect of noise 
delayed switching is clearly visible, quite similar to the 
one observed before for Josephson junctionsi^, and in 
general case of nonlinear systemsi^. Besides, here the 
SD behaves as tr ~ VT, see Ref J^. 

From Fig. 2 it is obvious that for 9=45° the reversal 
is faster and more stable than for 6'=0°,5° at all other 
equal conditions. Besides, the difference between MRT 
for the angles 9=0° and 9=5° is two times, while the 
difference of SD is about three times. For the cases 
9=5° and 0=45° the gain is even larger, more than 
one order for MRT and almost two orders for SD. This 
means that the reversal process principally depends on 
the precession of the magnetic dipole, and can not be 
described by a simple two-state model. This result gives 
the quantitative substantiation for the idea to use the 
tilted magnetic field to speed up the reversal processii, 
and also to use additional weak perpendicular magnetic 
field for the same purpose—, which actually leads to the 
tilt of the aggregate magnetic field. To roughly esti- 
mate the probability of nonswitching of a dipole by one 
pulse with the duration tp, the probability density of 
switching times can be considered as Gaussian with the 
mean r and SD a. Then the probability of nonswitching 
is: p=ierfc((tp - T)/V2a)). For T=300 K at the mini- 
mum of a we get p=10~'^^ even for 0=0°. However, for 
0=0°, 5° at the minimum of t we get unacceptably large 
probabilities 0.001 {p=10-^^ for 0=45°). 

In Fig. 3 the MRT and SD are presented versus the 
angle between anisotropy axis and external magnetic 
field for different values of damping a, tp=8 ns and 
T=300 K. For ^0, smaller values of damping lead to 
larger values of both r and a as it must (large values of 
T and a for —^0 and a=0.01 mean that in this range 
of parameters the noise- induced reversal occurs). For 
larger angles, however, the MRT nearly coincide, while 



the SD is smaller for smaller values of a. To understand 
why it is so, let us plot the trajectories of the magnetiza- 
tion for 0=0°, 5°, 45°, see the inset of Fig. 3 for a=0.1, 
tp=15.7 ns. One can see that the precession is largest for 
= 0°, for 0=5° the number of turns is smaller, and for 
0=45° the crossing of the boundary [0, My, Mz] occurs 
even without precession. This explains why the reversal 
in the latter case happens much faster than for 0=0° 
and has little dependence on a in the limit a<Cl, see 
Fig. 3. Since the length of the path is nearly the same 
for different a (the MRT nearly coincide) and the noise 
intensity is proportional to the damping, this obviously 
leads to smaller SD for smaller a. 

In Fig. 4 the MRT and SD versus angle are pre- 
sented for different values of the external magnetic field 
amplitude. It is seen that starting from values 20° — 30° 
and up to ~ 70° there are fiat minima of MRT and SD 
that corresponds to the range of angles where the fastest 
and the most reliable reversals are realized. For the tem- 
perature T = 300K the reversal occurs even up to the 
amplitude value Hq = 3kOe, which is two time smaller 
than the static coercive field He. 

In the present paper the effect of noise on the reversal 
of a magnetic dipole has been investigated on the basis 
of computer simulation of the Landau-Lifshits equation. 
It has been demonstrated that at the reversal by the 
pulse with sinusoidal shape, there exists the optimal du- 
ration, which minimizes the mean reversal time (MRT) 
and the standard deviation (SD, jitter). Also, both the 
MRT and the jitter significantly depend on the angle 
between the reversal magnetic field and the anisotropy 
axis. At the optimal angle the MRT can be decreased 
by 7 times for a=l and up to 2 orders of magnitude for 
a=0.01; the jitter can be decreased from 1 to 3 orders 
of magnitude (for a from 1 to 0.01) in comparison with 
the uniaxial symmetry case. For optimal angles the SD 
decreases with decrease of the damping. It has been 
demonstrated that fluctuations can not only decrease 
the reversal time, as it has been known before for the 
magnetic systems and is correct for small angles only, 
but it can also significantly increase the reversal time. 



* |alp@ ipm.sci-nnov.ru| 

^ A. Aharoni, Introduction to the theory of ferromagnetism, 

Oxford University Press, 2000. 
^ L. He, W.D. Doyle, and H. Fujiwara, IEEE Trans. Mag. 

30, 4086 (1994). 
^ A. Lyberatos, G. Ju, R. J. M. van de Veerdonk, and D. 

Weller, Journ. Appl. Phys., 91 2236 (2002). 

* D. Weller, and A. Moser, IEEE Trans, on Magn., 35, 4423 
(1999). 

^ H.N. Bertram, Wang Xiaobin, and V.L. Safonov, IEEE 

Trans, on Magn., 37, 1521 (2001). 
^ K.-Z. Gao, E.D. Boerner and H.N. Bertram, Appl. Phys. 

Lett., 81 4008 (2002). 

K.-Z. Gao and H.N. Bertram, Journ. Appl. Phys., 96 7753 
(2004). 

® M. Bauer, J. Fassbender, B. Hillebrands, and R. L. 

Stamps, Phys. Rev. B, 61, 3410 (2000). 
® G. Brown, M. A. Novotny, and P. A. Rikvold, Phys. Rev. 

B, 64 134422 (2001). 



Y. Uesaka, H. Endo, Y. Nakatani, N. Hayashi, and H. 
Fukushima, IEEE Trans, on Magn., 42, 1892 (2006). 
" A.N. Malakhov and A.L. Pankratov, Adv. Chem. Phys. 
121 357 (2002). 

A. L. Pankratov and B. Spagnolo, Phys. Rev. Lett. 93, 
177001 (2004); A.V. Gordeeva and A.L. Pankratov, Appl. 
Phys. Lett. 88, 022505 (2006). 

13 W.T. Coffey, D.S.F. Crothers, Yu.P. Kalmykov, and J.T. 

Waldron, Phys. Rev. B 51, 15947 (1995). 
" D.A. Garanin, Phys. Rev. E 54, 3250 (1996). 

Yu.P. Kalmykov, J.L. Dejardin, and W.T. Coffey, Phys. 

Rev. E 55, 2509 (1997). 

N.V. Agudov and A.N. Malakhov, Int. J. Bifurcation 
Chaos Appl. Sci. Eng. 5, 531 (1995); R.N. Mantegna and 

B. Spagnolo, Phys. Rev. Lett. 76, 563 (1996). 

Y. Y. Zou, J. P. Wang, C. H. Hee, and T. C. Chong, Appl. 
Phys. Lett., 82, 2473 (2003). 
1* R.L. Stamps and B. Hillebrands, Appl. Phys. Lett. 75, 
1143 (1999). 



